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Probing nonlinear mechanical effects through electronic currents: 
the case of a nanomechanical resonator acting as electronic transistor 

A. Nocera 1 ' 2 , C.A. Perroni 3 ' 4 , V. Marigliano Ramaglia 1,4 and V. Cataudella 3,4 
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We study a general model describing a self-detecting single electron transistor realized by a sus- 
pended carbon nanotube actuated by a nearby antenna. The main features of the device, recently 
observed in a number of experiments, are accurately reproduced. When the device is in a low 
current-carrying state, a peak in the current signals a mechanical resonance. On the contrary, a 
dip in the current is found in high current-carrying states. In the nonlinear vibration regime of 
the resonator, we are able to reproduce quantitatively the characteristic asymmetric shape of the 
current-frequency curves. We show that the nonlinear effects coming out at high values of the 
antenna amplitude are related to the effective nonlinear force induced by the electronic flow. The 
interplay between electronic and mechanical degrees of freedom is understood in terms of an unifying 
model including in an intrinsic way the nonlinear effects driven by the external probe. 



I. INTRODUCTION 



It has been recently shown that carbon nanotubes can 
act simultaneusly as single electron transistors^ (SET) 
and as nanoeletromechanical systems (NEMS/^ 3 . The 
idea is to use a single carbon nanotube placed between 
two metal contacts in a suspended configuration as a self- 
detecting SET. Due to the extreme properties of carbon 
nanotubes (ideal for NEMS applications, because they 
have a low mass and a high Young's modulus), the elec- 
tronic current flowing through the device results very sen- 
sitive to the dynamics of the nanotube itself. 

Although the main effort has been focused on detect- 
ing the quantum regime of mechanical resonators^"^, re- 
cently, G. A. Steele et affi and A. K. Huttel et a/P were 
able to fabricate a carbon nanotube electromechanical 
device working in the semiclassical regime (resonator fre- 
quencies in MHz range compared to an electronic hop- 
ping frequency of the order of tens of GHz) with an ex- 
tremely large quality factor (Q > 10 5 ). By measuring 
the variations of the electronic current flowing through 
the nanotube as function of the frequency of a nearby 
antenna actuating its motion, they were able to detect 
very well defined resonances corresponding to the bend- 
ing mode of the nanotube itself. The possibility of using 
currents to probe nanomechanical resonances is strongly 
related to the extremely large quality factor obtained. 
In particular, when the SET is in a low current-carrying 
state (far from electronic resonance), a peak in the cur- 
rent signals a mechanical resonance, while a dip is ob- 
served in a high current-carrying state (SET in elec- 
tronic resonance). Moreover, by adjusting the antenna 
power, the nanotube resonator can easily be tuned into 
the nonlinear vibration regime. For small antenna am- 
plitudes, the current-frequency curves are very well fitted 
by a Lorentzian, while a characteristic triangular shape 
or even hysteresis for large antenna amplitudes is ob- 
tained. Interestingly, the operating temperature can af- 



fect the nonlinearity and the quality factor of the res- 
onator in a non expected way: the nonlinear effects in 
current-frequency curves are washed out increasing the 
temperature^. A detailed analysis of a single frequency 
resonance dip has also shown that a broadening can be 
obtained increasing the source-drain voltage^. In Ref. 
[3 the authors were able to provide an explanation of 
some of the observed effects in terms of a model in which 
the gate voltage acquires an assigned time dependence. 
Within this phenomenological model, the back-action of 
the nanotube motion on detected current is understood 
neglecting completely the dynamics of the resonator and 
analyzing the problem directly at mechanical resonance 
conditions only in the limit of small external antenna 
amplitudes (linear response regime). 

In this paper, we do not consider a priori assump- 
tions on the resonator (nanotube) dynamics. Actually, 
in the device investigated in Ref. [7] and the chem- 
ical potentials of the leads, to which the nanotube is 
anchored, differ by the value of the applied transport 
voltage eVbias (where e is the electron charge), and thus 
the environment that the nanotube experiences cannot 
be considered at equilibrium since the voltages applied 
in the experiment are typically greater than the temper- 
ature, eVbias >> ksT (ks being Boltzmann constant). 
Therefore, in order to understand the behavior of the de- 
vice under such conditions, one needs to determine self- 
consistently the influence of electrical current and of the 
external antenna on the resonator dynamics, and vice 
versa, the influence of nonthermal nanotube vibrations 
driven by the antenna on the electronic currentffl As 
rigorously demonstrated by Mozyrsky et al\^, and re- 
obtained by us in a different way^l, in absence of the 
external antenna, the vibrational dynamics of the nan- 
otube can be described, employing a separation between 
slow vibrational and fast electronic time scales^ES by 
a Langevin equatiorp3ESl_ This equation is ruled by an 
effective force as well as a damping and diffusive terms 
stemming from the interaction of the resonator with the 



electronic bath consisting of both the nanotube itself 
(that can be described by a single electronic leveP^) and 
the out-of-equilibrium environment given by the macro- 
scopic leads. By including the external antenna effects 
through a forcing term in the Langevin equation, we pro- 
vide such a self-consistent description of vibrational and 
electronic dynamics as argued above. 

A theoretical treatment of the nanotube based device 
investigated in[7Jand[Hl including the external antenna ef- 
fects, has been already considered by G. Labadze and Ya. 
M. BlanteJ^l In Ref. [201 the authors use a Fokker-Plank 
equation for the resonator distribution probability based 



on master equation: 



This approach implicitely as- 



sumes that the energy scale of the applied voltages is 
much larger than the electronic tunneling energy scale 
for, eVbias » for. Moreover, tunneling is considered in 
the sequential regime and quantum effects in the elec- 
tronic dynamics such as cotunneling are disregarded^. 
We point out that these effects can be important to in- 
terpret the experimental results obtained by G. A. Steele 
et al™ and A. K. Huttel et alW since, as already em- 
phasized by G. Weick et al\^ 7 the temperature is much 
smaller than hT and the bias voltage effectively applied to 
the electronic level of the nanotube can be less than or of 
the same order of magnitude of the electronic tunneling 
energy, fc^T << eV^/J s < h~V^. Our approach, based on 
an adiabatic expansion of the time-dependent electronic 
Green function on the Kcldysh contour in the small pa- 
rameter wo/P^HIS takes into account from the beginning 
of all higher order terms in the tunneling matrix element 
between the leads and the nanotube (so sequential tun- 
neling and cotunneling regimes are described in a unified 
way). 

The inclusion of the external antenna effects in the ef- 
fective Langevin equation for the resonator has allowed 
us to explore also the nonlinear renspose regime in a 
nonperturbative way. We point out that, even describ- 
ing the vibrational dynamics of the nanotube as a sin- 
gle harmonic vibrational mode, we are able to reproduce 
all the main features observed in Ref. [7J and [5] includ- 
ing the nonlinear effects. Indeed, in the nonlinear vi- 
bration regime of the resonator, we are able to repro- 
duce quantitatively the characteristic asymmetric shape 
of the current-frequency curves. The observed nonlinear- 
ity is due to the intrinsic nonlinear terms of the effective 
force stemming from the strong interaction between elec- 
tron and resonator dynamic d 11 ' 12 * 17 * ^. These terms are 
completely neglected in Ref. [201 loosing any possibil- 
ity to describe renormalization frequency effects and to 
explore the nonlinear response regime for the resonator. 
We show that nonlinear effects can be highlighted only 
dynamically by applying large external antenna ampli- 
tudes. In this sense, our approach is different from other 
theoretical studies where nonlinear terms are discussed 
only statically-^ or are added from the beginning assum- 
ing t hat the resonator is characterized by anharmonic 
terms! 16 ' 25 ' 28 ' 

Within our approach, the experimental results ob- 



tained in Ref. [7J and [5] in the linear response regime are 
also reproduced, as well as the broadening of the mechan- 
ical resonance dip as function of the applied bias voltage. 
Furthermore, we are able to predict, in the limit of large 
bias, the onset of a fine double dip structure that could 
be experimentally observed. 

The paper is organized as follows: In Sec. II we present 
the model able to describe the electronic transistor con- 
sisting of the vibrating nanotube. The equation of mo- 
tion describing the resonator dynamics including the ex- 
ternal antenna effects is also discussed. In Sec. Ill we 
present numerical results. 



II. MODEL AND METHOD 

We describe the suspended carbon nanotube with a 
single impurity Holstein modePS that is able to catch 
the main physical ingredients of the experiments in Ref. 
[Jj and O For the sake of clarity, we here point out 
that the m odel us ed irFEl and in other papers in the 
literature^ * 23 * 25 '^, based on on a capacitive coupling of 
the nanotube (dot) to the gate electrode, is equivalent to 
a Holstein-like coupling between the occupation on the 
dot and the vibrational degree of freedom (see Appendix 
A). 

As suggested by Weick et al. [35] , in the small energy 
window of interest for a single dip feature, the electronic 
part of the device is modeled as a single electronic level 
coupled to the leads through standard tunneling terms. 

The electronic Hamiltonian is given by 



Hel = Vgl t {$ 'd + ^2(V k , a cl a d + h.c.) +y :/,.,/•;. ,/•/,-.,,. 



fc,« 



C fc,Q L 



eff 



(1) 

where nanotube's electronic level has energy V^/ e with 
creation (annihilation) operators d^df^*. The operators 
el a (ck, a ) create (annihilate) electrons with momentum 
k and energy Sk, a — Ek. a — Ha in the left (a = L) or right 
(a = R) free metallic leads, while the electronic tunneling 
between the molecular level and a state in the lead has 
amplitude Vh, a - The chemical potentials in the leads \ii, 
and /ifl are assumed to be biased by an external voltage 
e ^Uas = f-L — fJ-R- The coupling to the leads is described 
by the tunneling rate T a< k = 27rp a |V4 !Q ,| 2 /?i, where p a is 
the density of states in the lead a. We will suppose sym- 
metric coupling {Th,k = ^R,k) and a fiat density of states 
for the leads, considered as thermostats at finite temper- 
ature, wi thin t he wide-band approximation (r Q j. i— > T a , 

The Hamiltonian of the mechanical degree of freedom 



is given by H 



TTUOqX , characterized by the 



frequency wo and the effective mass m (k = mw 2 ). The 
interaction is provided by H int = \xN e p^2l, where A is 
the electron-oscillator coupling strength and N e i = <vd 
represents the electronic occupation on the nanotube (see 



also Appendix A). Definitely, the overall Hamiltonian is 



H = H e i+ H 



Hi 



(2) 



In this paper, we will measure lengths in units of xq = % 
and energies in units of ?i.r. 

The experimental values of the resonance frequencies 
of the vibrating nanotube (120-300 MHz range) sug- 
gest that the vibrational motion is very slow compared 
to the electronic tunneling rate on the nanotube itself 
(adiabatic limit): ujo/T << 1. Moreover, it was esti- 
mated in Ref. |25j that for the experiment in considera- 
tion the coupling energy describing electron-phonon in- 
teraction, E p = ^ ~ 5/j.eV, implying a strong coupling 
between the electronic and vibrational degrees of freedom 
(Ep/hujQ = 10). Summarizing, the regime of the relevant 
parameters is hoj a « E p (~ k B T) « eV^/J, < fiJ^. 

As discussed in Ref. [12] and [34], when eV^J a » hujQ 
and k B T > fiuo, the semi-classical treatment of the oscil- 
lator dynamics is well justifi ed. Wit hin a non-equilibrium 
adiabatic approximation^" ^ 13 * 18 * 19 *, the vibrational dy- 
namics of the nanotube can be described by a nonstan- 
dard Langevin equation controlled by a self-consistent ef- 
fective anharmonic force as well as by damping and dif- 
fusive terms depending explicitely on the resonator dis- 
placement x. The oscillator dynamics can be obtained 
by solving numerically the following equation 

mx + A(x)x = F(x) + £(t) + A ext cos(uj ext t), (3) 

where A ext and io ex t represent the amplitude and the 
external antenna frequency, respectively. Furthermore, 
in eq.©, the position-dependent force F(x), damping 
A(x), and the intensity of the noise D(x), (£(£)£(*')) = 
D(x)6(t — tf) (where £(£) is a standard white noise term), 
are related to the electroni c Greens functions on the 
Keldysh contour a d 11 ' 1 ^ ' 18 ' 19 ', 



F(x) = -kx + X{N el )(x) 

(Nei)(x)= ff^G<(u;,x) 
J 2m 



(4) 
(5) 



A(x) = A 



1 I d ^G<(u,x)d nu G>(uj,x), (6) 



D ( ,. )= A 2 / d ^G < (u,x)G > (u},x). 

2tt 



(7) 



The lesser G < and greater G > Greens functions at finite 
temperature are given by 



G < (u,x) 

G>(oj,x) 
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(8) 
, (9) 



where /l,.r(w) are the Fermi functions of the leads and 
we have assumed that KTl = TiTr = KT/2. At low tem- 
peratures, k B T << HT, the Fermi function f a (to) can be 



replaced by the step function <d(huj — fi a ) (a 
obtaining 

1 v^ ( f^a-V^-Xx 



L,R), 



(Nei)(x) = 2~ X! ( arctan ( 
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a=R 



The linear elastic force exerted on the oscillator is modi- 
fied by a relevant nonlinear correction term proportional 
to the electronic occupation eq.(|To|. The strength of 
the damping A{x) and diffusive D(x) terms result pro- 
portional to the adiabatic ratio ojq/T and therefore one 
can safely neglect their spatial dependence. Anyway, we 
point out that the diffusive term D(x) vanishes at equi- 
librium (bias voltage V^/J s = 0). Only on application of 
finite bias it becomes different from zero. In the regime 
of experimental parameters explored in Refs. [7J and |H1 
huj « E p (~ k B T) « eV b e /J s < hT, the effect of the 
electronic bath on the resonator dynamics can be de- 
scribed by an effective temperature proportional to the 
bias voltage k B T eff ~ e V£{{/d™ 

In eq.pl, we consider the parameter A ext expressed in 
terms of the natural force unit A = Lo^^/[2mE p ). As- 
suming a nanotube mass of m ~ 10 _23 fc<7, an oscilla- 
tion frequency of 120 MHz, A is of the order of 10~ 16 iV, 
while the effective spring constant is k = 10~ 6 N/m. By 
solving the second order stochastic differential equation 
(|3p2l, we are able to obtain the distribution probabilities 
of the displacement P(x). This allows us to calculate 
any system property as an average over the distribution 
probability P(x). In particular, in order to make contact 
with experimental results, we have calculated the aver- 
age electronic current (I) flowing through the nanotube 
as (I) = f_ dxI(x)P(x), where I(x) is the current at a 
particular resonator displacement. 



III. RESULTS 

One of the main results of Ref. [7J and |5] is the observa- 
tion of the electronic current changes at fixed gate voltage 
as function of the external antenna frequency and ampli- 
tude. Our model accurately reproduces the experimental 
results both in a low current-carrying state, where a peak 
in the current signals a mechanical resonance (FigjTa), 
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FIG. 1. (Color online) Panel(a): Normalized current change 
(AI/Iq) in a low current-carrying state (V e / t { = — 4XT) as 
function of the external frequency (uiext/uo) for different an- 
tenna amplitudes: A ex t 



line), 

A ex t 



10 solid thin line (red on- 
A ex t = 10 -4 ' 5 normal thickness line (green online), 
= 10" 3 - 5 thick line (blue online). Panel (b): AI/I 

ff 



against u> ex t/u>o in a high current-carrying state (Vg* te = E p 



for different antenna amplitudes: A e 



= 10" 



solid thin 



line (red online), A ext — 10~ 3 normal thickness line (green 
online), A ext = 10 -2 ' 5 thick line (blue online). Insets: solid 
line (black online) is a distribution P(x) out of mechanical res- 
onance. Dotted line (blue online) is a distribution obtained at 
mechanical resonance for the larger value of antenna ampli- 
tude considered in main plots of panels (a,b). Short-dashed 
line (magenta online) represents current as function of posi- 
tion I(x). In this plot eV b e / a { = 0.5KT, u /T = 0.004 and 

Ep/nr = k B T/hv = 0.04. 



and in a high current-carrying state, where a dip is ob- 
served (FigfT]}). 

These two behaviors (peak and dip) can be understood 
considering that the oscillator explores wider regions in 
configuration space when the external antenna amplitude 
increases (see the distribution probabilities P(x) shown 
in the insets of Fig. [I]) . When the electronic device is in a 
low current-carrying state (FigfT^i), P(x) is concentrated 
at x values far form the configurations where the device 
carries the maximum current (inset of FiglTk). By in- 
creasing the external antenna amplitude, the resonator is 
able to explore larger regions which carry more and more 
current obtaining a positive contribution in the normal- 
ized electronic current change AI/Iq with respect to the 
background value Iq. On the contrary, when the device is 
in a high current-carrying state, the distribution proba- 
bilities and the current are centered at the same position 
(inset of FigjTb). In this case, the effect of the external 
antenna is to give a negative contribution in the normal- 



ized electronic current change AI/Iq (FiglTh) since, the 
resonator explores regions of phase space which carry less 
and less current (inset of Fig. lib). 

The results shown in Fig[l] allow also to characterize 
the behavior of the resonator in the nonlinear regime. 
Interestingly, (see FigfT^i), increasing the amplitude of 
the external forcing, the shape of the current-frequency 
curves changes. For small antenna amplitudes, a char- 
acteristic Lorentzian shape is observed. This is expected 
for an harmonic oscillator driven by a periodic forcing 
in the absence of external noise. In fact, even if in me- 
chanical resonance, the oscillator explores only small re- 
gions around the stationary point and nonlinear correc- 
tions terms of the force F(x) (eq. Q) do not come into 
play. At mechanical resonance, only when the amplitude 
of the external antenna increases, the oscillator explores 
a larger region in phase space, where the nonlinear terms 
of the force acting on the oscillator cannot be neglected. 

Within our approach, for large antenna amplitudes and 
in the presence of noise, the current-frequency profiles 
assume the experimentally observed characteristic trian- 
gular shaped Furthermore, a softening is observed when 
the device is in a low current-carrying state (Fig. [1^), 
while an hardening in a high current-carrying state is 
obtained (Fig. ITp) . This nonlinear behavior can be un- 
derstood by analyzing the properties of the force F(x) 
(eq. El) around the stationary point. Softening and hard- 
ening behavior of the resonance frequency are usually 
related to the sign of the cubic nonlinear terrrpS. When 
the device is in a low current-carrying state, indeed, the 
sign of this term is positive, giving a net softening effect. 
In a high current-carrying state (V^ te = E p ) and for bias 
values sufficiently small, the sign of the cubic nonlinear 
term is negative providing an hardening. 

The temperature dependence of the current-frequency 
profiles, observed in the experiments, exhibits a non- 
trivial behavior that supports our model. As shown in 
Fig(2k, for very small temperatures, a triangular shape 
is found as discussed above. On the other hand, for suf- 
ficiently large temperatures, the current-frequency pro- 
file turns into Lorentzian shape, characteristic of the lin- 
ear response regime (Figl2p). This counterintuitive be- 
havior is determined by a significative reduction of the 
intrinsic nonlinear terms in the effective force F(x) on 
the resonator as function of the temperature. Actually, 
the correction term due the average electronic occupation 
(eq. (10)) tends to become independent of the resonator 
displacement x. Furthermore, the broadening of the 
current-frequency profiles as function of the temperature 
is produced not only by this effect but also by the growth 
of the intrinsic damping coefficient A(x) in the Langevin 
equation eq.J3]). Indeed, with increasing temperatures, 
A(x) increases where resonator dynamics occurs. The 
temperature dependence of the intrinsic damping A(x) 
is also responsible for the behavior of the resonator qual- 
ity factor Q, defined as 1$/ Alhalf-high- In the regime 
where the current-frequency profiles exhibit a Lorentzian 
shape, we find a power-law of the kind T~ 0A0 123 that 
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FIG. 2. (Color online) AI/Io against u) e xt/i^o in a low current- 
carrying state {Vgate ~ — 1.75Rr) when the resonator is driven 
by a strong external antenna amplitude A ex t = 10 -3 ' 5 : panel 
(a) k B T = 0.05M\ panel (b) k B T = 0.375fiT, panel (c) 
k B T = hF. In panel (c) a Lorentzian fit is also drown with 
Q = 10 4 (Q is defined in the main text). Inset of panel (a): 
intrinsic quality factor Q as function of temperature^. We 



use as energy unit hT = 125fieV. In this plot eV b 
w / r = 0.005 and E p /TiT = 0.05. 
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is quite close to that found in the experiment (T~ - 36 ). 
However, it is important to point out that this exponent 
depends on the gate voltage effectively applied to the 
device. 

Finally, we have analyzed all the traces of current vari- 
ations as function of the antenna frequency, obtained tun- 
ing the effective gate voltage. Going from a low to a high 
current-carrying state, the characteristic dip of the reso- 
nance frequency as function of the effective gate voltage is 
obtained in excellent agreement with experiments (curve 
(1) of FigJ3J). The observed renormalization of the reso- 
nance frequency can be related to strong variations of the 
electronic occupation eq.( 10 ) as function of the gate volt- 
age (see solid line in the inset of FigJ3J) . When the device 
is in a low current-carrying state (|V" g e / t { — E p \ > 1.57tr), 
the average electronic occupation is not sensitive to gate 
voltage variations. Instead, in a high current-carrying 
state (IF 



eff 

gate 



E p \ < 1.5HT), the electronic occupation 
shows a strong variation, providing the softening of the 
resonance frequencies. Increasing the bias voltage to val- 
ues closer to hT or larger (line (3) and (4) of Fig. 3J), 
one obtains a broadening of the resonance frequency dip, 
owing to a wider conduction window with respect to the 
broadening (~ hT) of the electronic energy level. Actu- 
ally, with increasing the bias, the electronic contribution 
to the effective spring constant increases, producing a 
nontrivial renormalization of the resonance frequency as 



FIG. 3. (Color online) Resonator frequency at resonance 
against effective gate voltage (shifted of E p ) for different bias 
voltages: eVbias = O.lhF curve (1), eVuas — Q.5KT curve 
(2), eV bl as = IMT curve (3), eVuas = 1.5ftr curve (4). 
Solid (red online) and short-dashed (blue online) portions 
of each curve indicate resonance frequency values with pos- 
itive and negative current change A/, respectively. Inset: 
electronic occupation at resonance frequency against effective 
gate voltage (shifted of E p ) for eVbias = 0.1KT (curve (1)) 
and eVbias = l.5hF (curve (4)). In this plot A ext = 10 -3 , 
ujo/r = 0.01 and E p /hT = k B T/hV = 0.1. 



function of the gate. We note that for eV^ ias — 1.57ir 
(line (4) of Fig. [3]), a fine structure represented by two 
very small dips appears. When the bias window, whose 
extension is proportional to eV^ , becomes larger than 
the broadening of the level, one could tune the electronic 
device into a region of conducting states where the varia- 
tions of the occupation are smaller than that obtained at 
the boundary of the conduction window itself (see dashed 
line in the inset of Fig 3l) . When the electronic device 
goes through states with different conducting character, 
the maximum renormalization of the resonance frequency 
occurs, providing two dips in the resonance frequency of 
the nanotube. This feature could be experimentally ob- 
served with a larger resolution in the applied gate voltage, 
when a very large bias is applied to the nanotube. 



IV. CONCLUSIONS 

In conclusion, we have studied a self-detecting SET re- 
alized by a suspended carbon nanotube including, in a 
non-perturbative way, the effect of the antenna driving 
the nanotube toward a nonlinear regime. All the qualita- 
tive features of the device, experimentally observed, are 
accurately reproduced clarifying the origin of the nonlin- 



ear effects. The nonlinear behavior is understood with- 
out adding by hand nonli near term s to the effective force 
exerted on the resonato j 7 1 16 1 25 1 28 1 , but stems out natu- 
rally from the nontrivial nonequilibrium time-dependent 
electronic occupation controlled by the coupling with the 
leads. We have shown that, increasing the temperature, 
the nonlinear effects in the current-frequency response 
are washed out as a result of the increase of the intrin- 
sic damping of the resonator and of the reduction of the 
intrinsic nonlinear terms of the effective self-consistent 
force. Within our approach, a broadening of the fre- 
quency dip as function of the bias voltage is reproduced, 
predicting in the limit of large bias a double dip struc- 
ture. 



SET dot and N gate = C ga teV ga te/e is the dimension- 
less electron number associated with a gate voltage V ga te 
which is coupled to the dot via a capacitance C ' ga te- 
In addition, a voltage Vuas is applied between source 
and drain which drives the tunneling of electrons across 



the SET. Here, C tot = C aa 



gate 



c, 



leads 7 



where C, 



leads 



is the sum of the capacitances resulting from the cou- 
pling to the leads. When the nanotube is allowed to 
vibrate, the gate capacitance C ga te assumes a spatial de- 
pendence C ga te(h(x)), where h is the distance between 
the nanotube and the gate electrode when the nanotube 
is displaced by a distance x from its equilibrium position 
(h(x) = ho + x). In the limit of small displacements, one 
can expand Hdot{x) around x — obtaining a Holstein- 
like linear correction term in x: 
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Appendix A 

Here, we briefly show that the model Hamiltonian 
for the vibrating nanotube encapsulated in the eq.pl) 
is equivalent to that u sed irpEl anc j j n other papers in 
the literature^HSMSEo] Indeed, as clearly discussed in 
Ref. [3H a single electronic transistor (SET) consists of a 
metallic dot (represented by a nanotube in the case con- 
sidered in Ref. [7] and |S]) with a large Coulomb-charging 
energy Ec = e 2 /2Ctot (Ctot is the total capacitance of 
the dot) coupled via tunnel junctions to both a source 
and a drain metal electrode. The Hamiltonian for the 
electronic and vibrational degrees of freedom of the dot 
is given by 



H int = AiV e iX, (A2) 

where A = —2(EcV g /e)(dC g /dx). In the small energy 
window of interest for a single dip feature investigated 
in Ref. [7] and we can neglect the weak gate and bias 
voltage dependence of A that is assumed constant! 21 * 32 * 33 !. 
Moreover, if we choose x = as the position where the 
Coulomb force for N e i = electrons in the island equals 
the elastic force, A can be then interpreted as the net 
force acting on the nanotube when one excess electron is 
populating the nanotube itself. The terms indepe ndent 
on N e i in the above expansion are usually neglecte d 32 * 33 !. 
One further assumes that the gate voltage is such that 
only charge states with N e i — and 1 are accessible. In 
this case, N^ = N e i and one can incorporate remaining 
constant terms (independent on x) in the above expan- 
sion in an effective gate voltage obtaining 
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dot 



v, 



e JJ e N el + XN el x + - 



gate 



kx 1 



(A3) 



H, 



dot 



Ec{N e 



Ngate 



r+^kx 2 , 



(Al) 



comprising a charging-energy term and the harmonic vi- 
brational energy. In eq.(Al), N e [ is the charge on the 



If we quantize the electronic and vibrational degrees of 
freedom, eq.( A2 ) gives the Holstein coupling discussed in 
the main text. 
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